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Abstract 

The  paper  outlines  procedures  for  extending  the  de  Casteljau,  de  Boor  and  Aitken  al- 
gorithms in  such  a way  as  to  allow  the  construction  on  a Riemannian  manifold  of  curves 
analogous  to  Bezier,  B-spline,  and  Lagrange  curves.  These  curves  lie  in  the  manifold  and 
respect  intrinsic  geometry. 


1 Introduction 

Given  a sequence  of  points  in  a Riemannian  manifold  M we  describe  methods  for  extend- 
ing the  de  Casteljau,  de  Boor,  and  Aitken  algorithms.  These  methods  allow  construction 
of  corresponding  interpolating  or  approximating  curves  that  lie  in  the  manifold  and  re- 
spect intrinsic  geometry.  In  the  case  that  the  manifold  is  a sphere,  opportunity  for 
applications  exist  in  the  domain  of  geological  and  geographical  mapping,  for  instance 
the  creation  of  topographical  contour  lines  or  isotherms,  and  in  the  field  of  video  pro- 
duction, where  it  is  desirable  to  have  smooth  camera  trajectories  interpolating  fixed 
camera  positions.  For  higher  dimensional  manifolds  there  are  applications  in  the  field 
of  data  analysis.  For  the  case  of  a sphere,  there  is  an  extensive  literature  dealing  with 
the  general  problem  of  data  fitting,  and  a superb  review  can  be  found  in  Fasshauer  and 
Schumaker  [2].  Shoemake  [7]  uses  properties  of  quaternion  arithmetic  to  describe  curves 
on  the  unit  quaternion  sphere,  and  Levesley  and  Ragozin  [4],  using  techniques  differ- 
ent from  those  presented  in  this  paper,  describe  methods  for  Lagrange  interpolation  in 
differentiable  manifolds. 

The  techniques  described  in  this  paper  come  from  the  simple  observation  that  in  the 
de  Casteljau,  de  Boor,  and  Aitken  algorithms  one  may  formally  substitute  appropriately 
parametrized  geodesic  arcs  for  straight  line  segments.  These  ideas  are  introduced  in  detail 
in  the  next  section  in  the  context  of  the  blossoming  paradigm,  [6]  and  [3].  Unfortunately 
many  of  the  useful  properties  of  blossoms  depend  on  the  affine  structure  of  Euclidean 
space  which  in  general  has  no  counter  part  in  a Riemannian  manifold.  In  particular, 
geodesic  blossoms  may  be  neither  symmetric  or  multi-affine,  and  in  general  they  do  not 
possess  uniqueness  characteristics  common  to  the  Euclidean  blossom. 

For  ari  arbitrary  Riemannian  manifold  [1]  or  indeed  an  arbitrary  differentiable  2- 
manifold  embedded  in  K3,  it  may  not  be  possible  to  construct  unique  shortest  geodesic 
arcs  between  two  points.  However,  if  the  manifold  is  compact  or  in  the  case  that  the  two 
points  lie  in  a sufficiently  small  neighborhood,  such  arcs  are  known  to  exist.  But  even 
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then,  there  appears  to  be  no  general  method  that  allows  explicit  construction.  So,  the 
task  of  constructing  geodesic  blossoms  becomes  a study  of  special  cases  in  which  specific 
methods  can  be  set  forth.  For  the  general  case,  a discrete  variational  method  can  be 
used  to  obtain  good  approximations. 

In  Section  3 a few  specific  examples  are  discussed.  The  case  in  which  the  manifold 
is  a sphere  is  given  special  attention.  There  we  introduce  a variation  which  allows  the 
discussion  of  Archimedian  curves  which  are  constructed  by  substituting  Archimedian 
spirals  for  geodesics.  This  variation  allows  the  natural  construction  of  curves  that  lie  off 
the  sphere.  Although  the  spherical  geodesic  blossoms  are  neither  symmetric  or  multi- 
affine,  a simple  reparametrization  of  geodesic  arcs  results  in  spherical  blossoms  that  have 
all  desirable  characteristics.  Section  3 also  contains  a brief  discussion  of  the  problem  of 
finding  geodesics  in  developable  surfaces  and  in  surfaces  of  revolution. 

2 Preliminaries 

Let  M be  a C°°  Riemannian  manifold.  There  is  the  following  theorem  that  guarantees 
the  existence  locally  of  geodesics. 

Theorem  2.1  If  M is  a Riemannian  manifold,  Xq  € M.  Then  there  exists  a neigh- 
borhood V of  £o  and  £ > 0 so  that  if  x £ V and  v is  a non-zero  tangent  vector  at  x 
and  Hvsll  < e,  then  there  is  a unique  C°°  geodesic  a : (—2, 2)  — *■  M defined  on  the  open 

interval  (—2,2)  such  that  a(0)  = x and  ( ^ ) = vx. 

\ at ) t=0 

For  compact  Riemannian  manifolds  there  is  the  Hopf-Rinow  theorem  that  tells  us 
that  points  can  be  connected  by  geodesic  arcs. 

Theorem  2.2  (Hopf  and  Rinow)  If  a connected  Riemannian  manifold  M is  compact, 
then  any  pair  of  points  x and  y may  be  joined  by  a geodesic  whose  length  corresponds  to 
the  distance  in  the  manifold  from  x to  y. 

We  also  need  the  notion  of  geodesic  convexity  and  the  result  of  J.  H.  C.  Whitehead 
that  geodesically  convex  neighborhoods  exist  for  all  x € M. 

Definition  2.3  Given  a subset  X of  M and  a point  xq  G X,  X is  star  shaped  with 
respect  to  the  point  Xq,  if  for  every  x £ X there  is  a unique  shortest  geodesic  connecting 
xq  with  x which  lies  in  X. 

Definition  2.4  A subset  X of  M is  geodesically  convex  if  it  is  star  shaped  with  respect 
to  each  of  its  points. 

Definition  2.5  Given  a subset  A of  a geodesically  convex  set  X the  geodesic  convex 
hull  of  A is  the  smallest  convex  set  which  contains  A. 

Theorem  2.6  (J.H.  C.  Whitehead)  LetV  be  an  open  subset  of  a Riemannian  manifold 
M and  let  x G M , then  there  is  a geodesically  convex  open  neighborhood  XJ  of  x such 
that  U C V. 

Let  M be  a Riemannian  manifold  and  let  A be  a geodesically  convex  subset  of  M. 
Given  points  Pi  in  M we  describe  extensions  of  the  de  Castlejau,  de  Boor,  and  Aitken 
algorithms. 
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2.1  Riemannian  Lagrange  curves 

Let  M be  a Riemannian  manifold,  and  let  A = {P0,  P1}  ■ ■ Pn}  be  a subset,  of  a 
geodesically  convex  subset  X.  Given  parameter  points,  to  < U < • • • < tn,  assume  that 
A is  contained  in  a sufficiently  small  neighborhood  in  which  specified  geodesics  exist. 
For  0 < i < n — 1,  define  jj  : [to,  t„]  — * X to  be  the  unique  geodesic  parametrized 
so  that  7 j(ti)  = Pi  and  = Pi+i-  For  1 < r < n and  0 < i < n - r define 

7 [ : [to,  tn]r  — » X so  that  7[(«i,«2,-  • ■,  ur-i , •)  is  the  unique  geodesic  parametrized  so 
that  7[(«x,u2)- • qUr-i,  U ) = 7[-1(«i,u2,- • qiir-i)  and  7[(ui,u2)  • • -.Ur-i,  U+r ) = 
7ir+i  (ui,U2,-  • qUr-i).  The  function  7q  : [t0,  tn]n  — > X is  called  the  geodesic  Aitken 
blossom  associated  with  the  points  P*  £ X,  0 < i < n and  the  parameter  points,  to  < 
tl  <•'••<  tn  • If  A : [to,  in]  — > [to,  tn]”  is  the  diagonal  map  defined  by  A(u)  = 
(u,  u,  • ■ ■ , u) , the  geodesic  Lagrange  curve  associated  with  X and  the  points  Pi  is  the 

V v ' 

function  Tg  = 7q  o A. 

Theorem  2.7  If  Tq  : [to,  tn]  —>  M is  the  geodesic  Lagrange  curve  associated  with  the 
points  Pi  £ M,  0 < i < n,  as  defined  above,  then  Tq  (t,)  = Pi. 

Proof:  Observe  that  for  1 < r < n and  0 < i < n - r,  7”  depends  for  its  definition 
only  on  the  points,  Pj,  where  i < j < i + r.  If  n = 1,  and  we  are  given  points,  Po  and 
Pi,  the  result  follows  from  the  definition  of  70  . Inductively  assume  it  is  true  for  k < n. 
For  k — n,  if  i = 0,  by  definition 

rj(to)  = 'To  (to,  to,-  • • , 1 0)  = 7o~1(t0,t0  ,•  • -,fq)  = • • • = 7o(<o)  = Po 

n n— 1 


and  likewise  if  i = n,  rg(t„)  = 7o  (tn,tn,--  • ,t„ ) = 7^  \tn,tn,--  -,tn)  = ■■■  = 7 £(<„)  = 

S v ' V " 

n n— 1 

P„.  For  i ^ 0 and  i ^ n,  observe  that  the  geodesics  used  in  the  construction  of 
7o_1  and  7”-1  may  be  restricted  respectively  to  the  intervals  [t0,fn_i]  and  [fi,tn]  so 
that  7q-1  becomes  the  geodesic  Aitken  blossom  associated  with  the  points  Po,Pi,'  ■ •, 
Pn_iand  the  parameter  points  to  < h < • ■ ■ < tn_  1,  and  7”-1  becomes  geodesic  Aitken 
blossom  associated  with  the  points  Pi,  P2,  • • •,  P«  and  the  parameter  points  t\  < £2  < 

• • • < tn.  By  the  deductive  assumption,  7q  _1(L,ti,  ••  •,£*)  = P = 7i  _1(* i,U,- • -,ti), 

' v ' 

n— 1 n— 1 

and  consequently  7 o(U,U,  • • -,U,  •)  is  the  geodesic  connecting  7o~1(f;,ft,  • ■ • ,U ) with 

n— 1 n— 1 

7 i~1(ti,ti,-  ■ and  is  thus  the  constant  function,  7o  (ti,U,-  • -,ti,u)  = Pi  for  all 

n— 1 n 

v £ [to,  tn].  Thus  in  particular,  7p  (tj,tj,--  -,tj)  = rg(f;)  = Pj.  □ 
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2.2  Riemannian  Bezier  curves 

Following  the  previous  format  we  introduce  a Riemannian  version  of  the  de  Casteljau 
algorithm.  Accordingly,  let  X he  a geodesically  convex  subset  of  a Riemannian  manifold 
M.  Let  A = {P0,Pi,  • ■ •,  Pn}  be  a subset  of  X.  Define  yf  : [0,  lj  — *•  X by  7? (it)  = 
Pi-  For  1 < r < n and  0 < i < n — r define  y\  : [0,  l]r  — > X to  be  the  unique 
geodesic  with  the  property  that  yf(u1,u2,-  • -.Wr-i,  0)  = 7[“1(ui,u2, • • • ,ur-i ) and 
7[(iii,  u2,  ■ • 1)  = 7[_j711(iii,  u2,  •• The  function  7q  : [0,  l]n  — ► X is  called 
the  geodesic  de  Casteljau  blossom  associated  with  the  set  A.  If  A : [0,  1]  — s-  [0,  1]"  is 
the  diagonal  map,  the  geodesic  Bezier  curve  associated  with  X and  the  set  A is  the 
function  Fq  = y^  ° A. 

2.3  Riemannian  B-Spline  curves 

Given  A = {P0,  Pi,  ■ ■ Pn}  contained  in  a geodesically  convex  subset  A of  a Rieman- 
nian manifold  M,  and  given  knots  U <t2  < ■ ■ ■ < t2n,  define  7?  : [fi,f2„]  — ► A by 
7 f(t)  = Pi,  for  0 < i < n.  For  1 < r < n and  r <i  <n,  define  y\  : [fi;  tj+7l+1_r]’’  — > X 
to  be  the  unique  geodesic  with  the  property  that  yf(ui,u2,  U ) = 7 T~* (v,i,u2,  •• 

■,ur-i)  and  7[(iii,u2)-  • -,ur_  1,  ti+n+ !_r)  = 7[_1(ui,u2)-  • -,ur_i).  The  function 
7 n '■  [tn,  tn+i}n  — > A is  called  the  geodesic  de  Boor  blossom  associated  the  set  A. 
If  A : \tn,  fn+i]  — > [tn,  tn+i]n  is  the  diagonal  map,  the  geodesic  B-Spline  curve  associ- 
ated with  A and  the  points  Pi  is  the  function  T”  = 7”  o A. 

We  have  the  following  results,  which  follow  from  the  fact  that  both  the  geodesic  de 
Casteljau  and  the  geodesic  de  Boor  blossoms  are  constructed  from  successive  geodesic 
combinations  beginning  with  the  set  A = {Po,  Pi,  • • •,  Pn}. 

Theorem  2.8  Given  A = {Po,Pi,---,  P„}  contained  in  a geodesically  convex  subset 
of  a Riemannian  manifold,  if  yj}  : [0,  1]”  ->  A is  the  geodesic  de  Casteljau  blossom 
of  A,  then  7o  ([0,  1]")  is  contained  in  the  geodesic  convex  hull  of  the  set  A. 

Theorem  2.9  Given  A = {Po,  Pi,  • • •,  Pn}  contained  in  a geodesically  convex  subset 
of  a Riemannian  manifold,  ifyft  : [ tn , f„+i]”  —>  A is  the  geodesic  de  Boor  blossom  of 
A relative  to  a knot  sequence  t\  < f2  < • • • < t2n,  then  yf; ([tn,  tn+ 1]")  is  contained  in 
the  geodesic  convex  hull  of  the  set  A. 

Since  each  of  the  three  blossoms  are  constructed  successively  from  C°°  geodesics,  it 
follows  that  the  blossoms  and  their  restrictions  to  the  diagonal  are  also  of  class  C°°. 

Theorem  2.10  The  geodesic  Lagrange,  Bezier,  B-spline  curves  are  of  class  C°°  as  are 
each  of  their  corresponding  blossoms. 

3 Examples 

The  impediments  to  implementation  of  these  ideas  depend  on  the  manifold  in  question. 
In  all  cases  it  is  necessary  that  the  points  Pj  should  lie  in  a region  in  which  it  is  possible 
to  construct  geodesic  arcs  between  points.  The  problem  then  reduces  to  that  of  finding 
methods  for  such  constructions.  Even  in  cases  for  which  this  is  possible,  there  is  the 
additional  problem  that  many  of  the  desirable  properties  associated  with  B-spline  or 
Bezier  curves  in  R3  may  have  no  direct  analogs.  Many  properties  such  as  the  ability 
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to  subdivide  a curve  depend  on  the  blossom  being  symmetric  or  multi-affine,  and  for 
the  generalizations  presented  here,  this  is  seldom  true.  For  the  case  of  an  orientable 
2-manifold  embedded  in  R3,  there  are  in  many  cases  good  solutions  to  the  problem  of 
finding  geodesics,  but  different  classes  of  surfaces  lead  to  different  solution.  In  this  section 
we  mention  a few.  In  the  case  that  the  manifold  M is  the  2-sphere  S 2 a preliminary 
version  of  our  results  is  reported  in  [5]. 


3.1  The  sphere 

In  the  case  that  M = S2,  a small  alteration  to  methods  presented  so  far  allows  the 
consideration  of  curves  that  lie  off  the  sphere.  Given  points  Band  Q that  lie  off  the 
sphere  consider  radial  projections  to  points  P and  Q and  let  7 : [a,  6]  — ► S2  be  a 
geodesic  with  the  property  that  7(a)  = P and  7 (b)  = Q.  The  curve  7 : [a,  b]  — > R3 
defined  by 

7W  = (H-Hpll  + H-|IQI|)-7W 

is  called  the  Archimedian  spiral  connecting  the  points  P and  Q.  To  explicitly  describe 
the  curve  7 , set  P = uj,  Q = v2  and  for  simplicity  consider  the  parameter  interval  [a,  b ] 
to  be  the  unit  interval  [0, 1],  For  < •,  • > the  standard  inner  product  on  R3  set 

V3  = (<  VUV2  >Vi-  u2)/(||<  Vi,V2  >vi  - v2||) 

so  that  V3  is  orthogonal  to  Vi  and  in  the  plane  containing  v\  and  v2.  Letting  6 =<  v\,v2  > 
denote  the  angle  between  17  and  v2,  the  geodesic  7 connecting  v\  with  v2  is  defined  by 


7 (t)  = cos  (td)vi  + sin(t,9)v3 

sin(t$)  < vi,  v2  > 


cos  {t9)  + 


|<  Vi,V2  > Vi  - v2\\ 


vi  - 


sin(t0) 


|<  Vi,V2  > Vi  - v2\ 


;v2- 


The  corresponding  Archimedian  Lagrange,  Bezier  and  B-spline  curves  may  now  be  con- 
structed with  the  general  algorithms  of  Section  2. 

One  of  the  difficulties  that  arise  with  Archimedian  curves  is  that  geodesic  blossoms  are 
not  necessarily  symmetric  or  multi-affine.  It  is  even  not  clear  what,  these  concepts  might 
mean  in  a geodesic  context.  Consequently,  certain  results  that  hold  for  normal  Bezier  or 
B-spline  curves  that  depend  on  these  properties  are  no  longer  valid.  In  particular  analogs 
of  the  subdivision  algorithms  that  allow  one  to  determine  control  points  of  a portion  of  a 
given  Bezier  or  B-spline  are  not  valid.  However,  it  can  be  shown  that  a simple  non-linear 
change  in  the  parametrization  of  the  geodesic  arcs,  makes  it  possible  to  recapture  most 
of  what  is  needed. 

Definition  3.1  Given  two  points  A and  B on  the  sphere.  Let  C be  the  smaller  arc  of 
the  spherical  geodesic  joining  A with  B.  The  bary centric  param.etriza.tion  of  C on  the 
parameter  interval  [a,  6]  is  the  function  a : [a,  6]  — > C defined  by 


a{t)  = q(x(t)), 


In  the  following  we  prove  a spherical  version  of  the  Menelaus  theorem. 
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Theorem  3.2  Given  3 points  Pq,  P\,P2  on  Si 2  let  7 : [0, 1]  x [0, 1]  — > R3  be  the  geodesic 
de  Casteljau  blossom  in  which  all  geodesic  arcs  are  given  the  barycentric  parametrization. 
Then  j(s,t)  = 7(t,  s). 

Proof:  Observe  that  an  elementary  geometric  argument  tells  us  that: 

/y(s,t)  =7o(M)  = ?((! -i)7o(s)  + <7Ks)) 

= q(  (l-<)[(l-s)P0  + sPi]+f[(l-s)Pi+sP2]  ) 

and 


7(*,s)  = 7o(M)  = g((l-s)7o(0  + «7iW) 

= <?(  (1  - s)[(l  - t)P0  + tPi]  + s[(  1 - t)Px  + tP2)  ). 

And  the  result  follows  from  the  afhne  properties  of  K3. 


□ 


As  an  immediate  consequence  we  have 

Theorem  3.3  Given  points  Po,  Pi,  Pn  on  S2,  the  associated  de  Casteljau  blossom, 
in  which  geodesic  arcs  are  given  barycentric  parametrization,  is  symmetric. 

The  conventional  blossoming  description  of  subdivision  can  now  be  employed.  From 
the  blossom  construction  we  can  conclude  that  7J (0, 0,  - - •, 0,  1, 1,  ■ • •,  1)  = P,.  In  par- 

V*  v- 

i 

ticular,  it  follows  that,  for  0 < u < 1,  the  points  Q,  = 7o(0,0,  ■ • 0,  it,  it,  • ■ -,it) 

i 

describe  a geodesic  de  Casteljau  blossom  which  is  parametrized  to  the  interval  [0,tt] 
and  which,  because  of  the  uniqueness  of  geodesic  arcs,  equals  the  restriction  of  7#  to 
[0,u]n.  Likewise,  for  the  interval  [u,  1]  the  points  Ri  = 7 q(u,u,  •••,«,  1, 1,  • • •,  1)  de- 

i 

termine  a geodesic  de  Casteljau  blossom  which  is  parametrized  to  the  interval  [it,  1] 
and  which  equals  the  restriction  of  7J  to  [u,  1]”  Therefore,  if  g : [0, 1]  — *■  S2  is  the 
geodesic  Bezier  curve  determined  by  Pq.  Pi,  - ■ • , Pn  and  if  g = Jq  o A,  it  follows  that, 
s|[o,u]  = 1 7o(Mr  • ‘,t,  u, «,•••, it)  and  £t|[liU]  : 1 1->  7g(it,n,  • • -,it,  t,  t,---,t),  for 

i i 

0 < n < 1. 

More  generally  and  along  the  lines  of  the  proof  above,  we  have  the  following  theorem 
which  allows  all  familiar  properties  of  both  Bezier  and  B-spline  curves  which  have 
descriptions  in  terms  of  their  corresponding  blossoms  to  carry  over  to  the  spherical 
case. 

Theorem  3.4  Let  f : [0, 1]”  — > R3  be  the  Euclidean  blossom  generated  by  the  de 
Casteljau  algorithm  using  points  Pi  £ S2, 0 < i < n.  Then  7o  — 9 0 /• 

3.2  Other  surfaces 


We  briefly  discuss  two  examples  in  which  explicit  descriptions  of  geodesics  between  points 
are  possible. 

A developable  surface  S [4],  described  as  the  image  of  a function  / : i7  — > R3  for 
U an  open  subset  of  R2,  possess  the  characteristic,  among  others,  that  distances  are 
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preserved  by  the  function  /.  Therefore,  a geodesic  in  the  surface  f(U)  may  be  considered 
as  the  image  of  a straight  line  in  the  plane.  If  Pq,Pi,-  • yP„  are  points  in  S , let  Q,  = 
f~1(Pi),  0 < i < n.  If  C C U is  the  Lagrange,  Bezier,  or  B-spline  curve  obtained 
from  the  standard  Euclidean  versions  of  the  algorithms,  then  it  follows  that  f(C)  is  the 
corresponding  geodesic  curve  in  S that  would  have  been  obtairied  using  geodesic  versions 
of  the  algorithms  that  we  have  described. 

For  surfaces  of  revolution  the  description  of  geodesics  between  two  points  is  rather 
more  involved.  Let  C be  a curve  in  the  yz- plane  described  implicitly  by 

f f{y)  = z 

\ x = 0 ’ 

for  ( y , z)  belonging  to  some  open  set  U contained  in  the  upper  half  of  the  7/2-plane.  The 
surface  S obtained  by  rotating  C about  the  z-axis  may  be  expressed  as  g~1( 0)  where 
g : R X U — ► R is  defined  by  g(x,y,z)  = /(1/x2  + y 2)  -2  = 0.  In  polar  coordinates 
letting  u = \] x 2 + y2,  we  express  S in  the  form 

{x  = u cos  9 
y = u sin  9 . 
z = f(u) 


Let  P = (ui  cos 0i,  tii  sin /(tii))  and  Q = (772  cos #2,7*2 sin #2, /(«2))  be  two  points 
on  S.  Then  it  may  be  shown  that  the  geodesic  connecting  P with  Q is  the  function 
a : [711,7x2]  —>  S such  that  a(u)  = (u  cos 9 (v),  sin  9 (u),  f(u)),  where  for  fixed  txq, 


1 + (/W 


dt  + c! 


and  constants  c and  c'  satisfy  the  following  equations: 


9i  — 9\  = 


r ;i  + (/w 
L V - 1,2 


du 


'1  + (/M>2 


tv7i4 


du. 


For  complete  details  see  [6]. 


4 Conclusion  and  future  research 

We  have  outlined  a procedure  by  which  conventional  computer  aided  design  constructions 
may  be  extended  to  arbitrary  Riemannian  manifolds.  In  practice,  there  are  difficulties. 
In  a given  manifold  points  to  be  interpolated  or  approximated  must  lie  in  a region  in 
which  it  is  possible  to  construct  necessary  geodesic  arcs.  Supposing  this  the  case,  one 
then  needs  to  find  explicit  descriptions  of  the  geodesics.  And  then  there  is  the  question 
of  the  additional  characteristics  which  the  curves  might  possess.  The  paper  raises  more 
questions  than  it  answers.  In  the  case  of  a sphere,  good  results  are  obtained,  and  it 
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is  also  possible  to  add  variation  that  allows  consideration  of  curves  off  the  sphere  but 
which  project  radially  to  geodesic  Lagrange,  Bezier,  or  B-spline  curves.  It  is  also  shown, 
in  the  spherical  case,  that  a change  parametrization  of  geodesics  results  in  blossoms  that 
retain  the  desirable  characteristics  associated  with  Euclidean  blossoms.  For  surfaces  of 
revolution  and  developable  surfaces,  we  know  that  geodesics  can  be  found  between  points 
so  the  geodesic  blossom  constructions  will  always  exist.  It  is  however  unlikely  that  these 
blossoms  will  be  either  symmetric  or  multi-affine;  these  characteristics  depend  on  the 
affine  structure  of  R3.  Thus,  in  the  case  of  a general  Riemannian  manifold,  although  the 
constructions  may  be  valid,  it  is  not  clear  that  we  will  be  able  to  employ  fundamental 
operations  such  as  subdivision  which  depend  on  the  symmetry  of  the  blossom.  We  have 
outlined  three  different  methods  of  blossom  construction,  one  for  each  of  the  algorithms 
considered.  In  the  Euclidean  case,  we  know  that  there  is  a unique  symmetric,  multi- 
affine  polynomial  that  restricts  to  a given  polynomial  on  the  diagonal.  This  may  not  be 
true  in  our  more  general  setting. 
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